Jastrow theory of the Mott transition in bosonic Hubbard models 
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We show that the Mott transition occurring in bosonic Hubbard models can be successfully 
described by a simple variational wave function that contains all important long-wavelength corre- 
lations. Within this approach, a smooth metal-insulator transition is made possible by means of a 
long-range Jastrow correlation term that binds in real space density fluctuations. We find that the 
Mott transition has similar properties in two and three dimensions but differs in one dimension. We 
argue that our description of the Mott transition in terms of a binding-unbinding transition is of 
general validity and could also be applied to realistic electronic systems. 

PACS numbers: 71.10.Hf, 71.27.+a, 71.30.+h 
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Stimulated by the discovery of many strongly- 
correlated materials which, on the verge of becoming 
Mott insulators, display interesting and unusual prop- 
erties, a huge theoretical effort has been devoted in the 
last decades to clarify the interaction-driven Mott metal- 
insulator transition (MIT). 1] In spite of that, a full 
comprehension of this phenomenon is still lacking, even 
though, in the limit of infinite-coordination lattices, the 
whole dynamical behavior across the MIT can be un- 
covered thanks to Dynamical Mean-Field Theory. Q As 
a matter of fact, the Mott phenomenon is not specific of 
fermions but also occurs in bosonic systems, 0] that have 
recently become popular in the context of optical lattices, 
where a MIT can be actually realized experimentally. |4] 

The prototypical Hamiltonian to describe a MIT both 
for fermions and bosons is the Hubbard model 



H 



+ H.c.)+^J2 n >( n *- 1 )> W 



where b irT (bi a ) creates (annihilates) a particle at site i 
with spin a = —S,...,S, being S half-odd-integer for 
fermions and integer for bosons and rn = b\ a b ia the 
local density operator. The Hubbard model (JTJ) at com- 
mensurate densities generally shows two different phases: 
If U < U c particles are delocalized, which implies a 
metallic behavior for fermions (unless a Stoner instabil- 
ity leads to magnetically ordered phases) and superflu- 
idity for bosons; instead, when U > U c the model de- 
scribes a Mott insulator where coherent motion is sup- 
pressed. Presumably, for any S ^ 0, the Mott insulat- 
ing phase is accompanied by translational and eventually 
spin-rotational symmetry breaking. However, the latter 
is merely a consequence of the Mott phenomenon and 
should not be identified as the driving mechanism lead- 
ing to the insulating behavior, which arises from a strong 
suppression of charge fluctuations. Indeed, a MIT docs 
occur also for S = bosons, in which case no symmetry 
breaking is expected within the Mott insulator. In this 
work we focus on the bosonic Hubbard model by describ- 



ing the Mott phase with a strongly correlated variational 
wave function. An advantage of considering bosons is 
that we can directly compare the variational outcome 
with numerically exact results obtained by Green's func- 
tion Monte Carlo (GFMC). @, @] Therefore, by means 
of this comparison, we can establish the key ingredients 
that must be included in the variational state for a faith- 
ful representation of a genuine Mott insulating state. 

In spite of the fact that the variational approach is 
a simple and well established technique, its application 
to the MIT turns out to be extremely difficult. For in- 
stance, the celebrated Gutzwiller wave function \^g) = 
EL li n i)\^o) i where |$o) is the non- interacting ground 
state and "f(rii) is an operator which progressively sup- 
presses expensive occupancies, is not appropriate to de- 
scribe the MIT in both fermionic and bosonic cases. 
Indeed, the only way to produce an insulating wave func- 
tion corresponds to project out completely on-site oc- 
cupancies different from the average one. This wave 
function, with no charge fluctuations, is clearly a very 
poor description of a realistic Mott insulator. In fact, 
for the fermionic Hubbard model, the optimal |^g) is 
never insulating, except at U = oo, even if the variational 
wave function is improved by adding short-range density- 
density correlations. 0, 0] In the case of S = bosons, 
an insulating \^g) can be stabilized at finite U, § but, 
as we mentioned, the insulator obtained in this way gives 
an incorrect description of the actual ground state. 

A step forward has been recently accomplished in one 
dimension, where it has been shown that a Gutzwiller 
wave function supplemented by a long-range Jastrow fac- 
tor offers a very accurate description of a Mott insulator. 
However, a systematic analysis of this variational ansatz 
in higher dimensions is still lacking, while it would be 
highly desirable in view of more realistic applications. 
In this letter, we apply this variational approach to the 
5 = bosonic Hubbard model (p} with nearest-neighbor 
hopping t/2 in a one-dimensional chain (ID), a two- 
dimensional (2D) square lattice and a three-dimensional 
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FIG. 1: Variational results for the Jastrow potential v q mul- 
tiplied by q 2 in ID and 2D and by |g| 3 in 3D for increasing 
values of U/t (from bottom to top). Upper panel: ID case for 
60 and 100 sites. Middle panel: 2D case for 20 x 20, 26 x 26, 
and 30 x 30 clusters (along the (1,0) direction). Lower panel: 
3D case for 8 x 8 x 8, 10 x 10 x 10, and 12 x 12 x 12 clusters 
(along the (1,0,0) direction). 
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FIG. 2: Variational values of the condensate fraction n /L in 
ID (upper panel), 2D (middle panel), and 3D (lower panel). 



(3D) cubic lattice with L sites. Specifically, we consider 
the following ansatz for the variational wave function 



l*> = CX P {^\Y1 Vi 



QMB^ti i $ °>> ( 2 ) 



where |$o) is the non- interacting fully-condensed wave 
function, i.e. |$o) = (6fc =0 ) JV |0), being b\ the creation 
operator at momentum k and N the number of particles. 
In the following, we will consider N = L. The compo- 
nents of the Jastrow potential, v% j = u(|i2j — Rj\), are 
independently optimized by minimizing the variational 



energy, [ll[ and we will denote by n q and v q the Fourier 



transforms of the local density m and of Vij, respec- 
tively. Finally, Qmb is a variational parameter related to 
the many-body operator^ = hiYl s (l-d i+ s)+diYl s (l- 
hi+s), where hi = 1 {d% = 1) if the site i is empty (doubly 
occupied) and otherwise, and S is the vector which con- 
nects nearest-neighbor sites; [l2[ in other words, J2i & 
counts the number of isolated empty and doubly occu- 
pied sites. This term is kept mainly to improve the vari- 
ational accuracy (in 2D and 3D) but does not introduce 
important correlation effects, that are instead contained 
only in the long-range behavior of the two-body Jastrow 
potential Vij. Remarkably, it turns out that our wave 
function © is quite accurate in all cases that we consid- 
ered, even across the MIT; more details will be reported 
elsewhere. 
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Let us start by discussing the relevance of the Jastrow 
factor for the low-energy properties and some expected 
asymptotic behaviors of v q . In the gapless (superfmid) 
phase, a long-range Jastrow potential is surely needed to 
restore the correct small- g behavior of the static density 
structure factor, i.e., N q = (*|n_ g n g |*)/(*|*) ~ |g|. 
Indeed, since at least in the weak-coupling regime, the 
expression 



N„ 



(3) 



holds with 7 = 2, 14] where the non-interacting N° = 
(&o\n- q n q \&o) / ($>o\$o) ~ const, it follows that v q ~ 
l/|g|. The same asymptotic behavior was obtained in 
Refs. QI EE 

by full optimization of v q in metal- 
lic fermionic models both in ID and 2D, where both 
Moreover, it was shown that, in 



N„ 



\q\ and N° 



the insulating phase, a more singular v q ~ 1/q 2 at small 
q is required to recover the appropriate N q ~ q 2 insu- 
lating behavior, consequence of exponentially decaying 
correlation functions. This similarity between ID and 
2D fermionic models was suggestive that v q ~ 1/q 2 is 
sufficient to induce an insulating behavior in any dimen- 
sions as well as that the expression (j3|) remains asymptot- 
ically valid for \q\ — » 0, even inside the insulating regime 
at strong coupling. However, one easily realizes that, 
were this conclusion correct, the variational wave func- 
tion could not describe any bosonic insulator in 3D, 
since v q ~ 1/q 2 is not sufficient to empty the conden- 
sate fraction n = (*|4 =0 b k=0 \^) / {^\^}- [13] Instead, 
we will show that our Jastrow wave function can give a 
consistent description of the Mott phase in any dimen- 
sion thanks to an even more diverging v q , implying that 
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the formula ([3]) can be violated in the real 3D world. 

In Fig. [T] we present the optimized Jastrow potential 
V q . For all dimensions, the MIT can be clearly detected 
from the sudden change in the small-g behavior of v q . 
On the one hand, the gapless superfluid phase is always 
described by v q ~ with a increasing with U. On 

the other hand, the gapped insulator has a much more 
diverging v q . In ID we recover the v q ~ 1/q 2 behavior, 
like in the fermionic case. [l(| In 2D, the leading behav- 
ior of the Jastrow potential across the transition is less 
clearcut than in ID. Indeed, we cannot establish whether, 
on the insulating side, v q ~ fau/q 2 with P2D large but 
finite, or possible logarithmic corrections have to be con- 
sidered, i.e., v q ~ ln(l/|g|)/<7 2 . The first possibility is 
particularly appealing since, in this case, the insulating 
phase can be interpreted in terms of the confined phase 
of the 2D classical Coulomb gas. H, 18| Notice that the 
optimized v q , that also contains subleading corrections 
to the 1/q 2 behavior, can modify the critical properties 
of the classical Coulomb gas model. Nevertheless, the 
essential point is that, within this approach, the MIT 
can be still interpreted in terms of a binding-unbinding 
transition among charged particles (empty and doubly 
occupied sites). Finally, in 3D an even more diverging 
v q is stabilized in the insulating regime, i.e., v q ~ l/|g| 3 . 
Therefore, in all these cases the Jastrow potential is suf- 
ficient to destroy the condensate (see Fig. [2]). [l^ 

In order to verify the validity of our approach, let us 
move to discuss the variational results for the density 
structure factor N q in comparison with the exact ones 
obtained by GFMC. For small q's we can generally write 
N q = Ji\q\ + 72 <Z 2 + 0(g 3 ). In analogy with spin sys- 
tems, we may assume that 71 = v c x, being v c and x 
the charge velocity and the compressibility, respectively. 
At the variational level 71 and 72 depend crucially upon 
the Jastrow parameters. Indeed, we do find that, in the 
superfluid phase, 71 ^ while, in the Mott insulator, 

71 = and 72 ^ 0, signaling that this state is incompress- 
ible (see Fig. [3]). Moreover, in ID we have evidence for a 
jump (from a finite value to zero) in 71 across the MIT, 
especially because its value does not change much passing 
from U = U c /2 and U <U C (e.g., it changes from 0.4 to 
0.2). These variational results are confirmed by GFMC 
and are consistent with the finite jump of the compress- 
ibility across the Kosterlitz-Thouless transition, expected 
in ID. Q Our numerical results seem also to indicate that 

72 diverges as the MIT is approached from the insulating 
side. In the variational calculation, this behavior follows 
from a v q ~ f3uj/q 2 in the insulating phase with ftm — » 
at the MIT. In conclusion we find that the ID MIT can 
be located at U c /t ~ 2.45 in the variational calculations, 
whereas the GFMC gives U c /t ~ 2.2 (in close agreement 
with previous calculations of Ref . [2(], [2l[ ) , showing that 
the variational wave function (|2|) is not only qualitatively 
but also quantitatively correct. 

The density structure factor N q displays quite distinct 
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FIG. 3: Density structure factor N q divided by \q\ calculated 
with variational Monte Carlo (left panels) and GFMC (right 
panels) in ID (upper panels) and 2D (lower panels) In ID, 
L — 60 and U/t = 1.6, 1.8, 2, 2.2, 2.4, 2.5, 3, and 4. In 2D, 
L = 20 x 20 and U/t = 10, 10.2, 10.4, 10.6, and 10.8 for the 
variational calculation, and L — 256 and U/t = 8, 8.2, 8.4, 
8.6, and 8.8 in the GFMC calculation. All cases are shown 
from top to bottom for increasing values of U /t. 
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FIG. 4: Upper panel: Density structure factor iV ? divided 
by q 2 calculated by the variational Monte Carlo for 3D and 
U/t = 20. Lower panels: N q for non-optimized wave functions 
with v q ~ /33o/|g| 3 for two values of P3D and the same sizes 
as above. 
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long-wavelength behaviors for weak and strong interac- 
tions also in 2D, see Fig. [3] In the variational calcu- 
lations, for U/t < 10.3 the structure factor goes like 
N q ~ j t \q\, while for U/t > 10.3 we get N q ~ 7 2 g 2 . The 
critical value of the on-site interaction is slightly different 
from the GFMC one, which we find to be U c /t ~ 8.5, in 
agreement with Ref. 22j . In spite of the different values 



of U c , the qualitative behavior across the MIT is similar 
both in the variational and in the GFMC calculations. 
It should be emphasized that in 2D the value of 71 close 
to the MIT is one order of magnitude smaller than its 
value at U/U c ~ 1/2, a behavior qualitatively different 
from ID, suggesting that 71 vanishes upon approaching 
the MIT from the superfluid side and that 72 is smooth 
and constant across the MIT in 2D. 

More interesting is the 3D case. Here, the GFMC is 
severely limited by small sizes and, therefore, we will just 
discuss the variational results. The change in the lead- 
ing behavior of the Jastrow parameters v q allows us to 
locate the transition around U c /t ~ 18, which is very 
close to the critical value extracted from experiments on 
optical lattices. Q The optimal Jastrow potential be- 
haves as usual as v q ~ a/\q\ in the superfluid phase, 
but turns into v q ~ /?3.d/|<z| 3 in the Mott insulator (see 
Fig. Q]) • This behavior would imply, if Eq. ((3J were valid, 
a charge structure factor N q ~ \q\ 3 . By contrast, we 
do find that N q ~ q 2 , as expected in an insulator, see 
Fig. HI So, we arrive at the very surprising and unex- 
pected result that Eq. ([3]) does not hold, not even asymp- 
totically for \q\ — ► 0. In order to prove more firmly that 
a v q ~ /?3Z? /|f?| 3 can indeed lead to N q ~ q 2 , we have 
calculated N q with a non-optimized wave function of the 
form ^ with v q ~ @3D/\q\ 3 , for different values of /?3_d. 
As shown in Fig. [4] for small /33d's N q ~ |q| 3 , implying 
that Eq. ([3]) is qualitatively correct. However, above a 
critical (3^ D , the behavior turns into N q ~ q 2 , signaling a 
remarkable breakdown of Eq. ([3]). The optimal value of 
@3D that we get variationally at the MIT is larger than 
/?3£), confirming our variational finding N q ~ q 2 . We no- 
tice that the change of behavior as a function of (3%d is 
consistent with the binding-unbinding phase transition 
recently uncovered in a classical 3D gas with potential 
V d (q) ~ l/\q\ 3 . M 

In conclusion we have demonstrated that a long-range 
Jastrow potential does allow for a faithful variational de- 
scription of a Mott transition in the bosonic Hubbard 
model, in spite of the fact that the uncorrelated wave 
function onto which the Jastrow factor is applied has full 
Bose-condensation. An interesting outcome of our analy- 
sis is that, in 3D, the Mott insulator is characterized by a 
very singular Jastrow potential, v q ~ l/|g| 3 , that is able 
to empty the condensate, yet leading to a well behaved 
charge structure factor, N q ~ q 2 . This result contra- 
dicts the naive expectation, N q ~ l/v q , based on the 
weak-coupling formula ^ . This breakdown of the weak- 
coupling approach is the necessary condition for our wave 



function to work in 3D and represents a highly nontriv- 
ial consistency check of our non-perturbative variational 
theory of the Mott phase. We argue that this variational 
theory will hold also in electronic models. In particu- 
lar, once the square of the ground state wave function is 
interpreted as a classical partition function, the metal- 
insulator transition can be induced in any dimension by 
a singular interaction between charge fluctuations. Re- 
markably, in D > 1 this interaction remains always loga- 
rithmic, suggesting an unconventional binding-unbinding 
description of the metal-insulator transition. In analogy 
with the bosonic example we have analyzed, we should 
expect that a singular Jastrow potential, v q ~ ^/\q\ e with 
= 3, might be necessary to describe the 3D Mott tran- 
sition in fermionic models, too, all the more reason when 
realistic Coulomb interaction is taken into account. 
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